A method of alternating characteristics with application to advection-dominated environmental systems     

Katerina Georgiou  John Harte  Ali Mesbah  William J. Riley 《Computational Geosciences》2018,22(3):851-865

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.  相似文献   

Seasonal variability of oxygen and hydrogen isotopes in a wetland system of the Yunnan-Guizhou Plateau,southwest China: a quantitative assessment of groundwater inflow fluxes     

Xingxing Cao  Pan Wu  Shaoqi Zhou  Han Tu  Shui Zhang 《Hydrogeology Journal》2018,26(1):215-231

The Caohai Wetland serves as an important ecosystem on the Yunnan–Guizhou Plateau and as a nationally important nature reserve for migratory birds in China. In this study, surface water, groundwater and wetland water were collected for the measurement of environmental isotopes to reveal the seasonal variability of oxygen and hydrogen isotopes (δ¹⁸O, δD), sources of water, and groundwater inflow fluxes. Results showed that surface water and groundwater are of meteoric origin. The isotopes in samples of wetland water were well mixed vertically in seasons of both high-flow (September) and low-flow (April); however, marked seasonal and spatial variations were observed. During the high-flow season, the isotopic composition in surface wetland water varied from ?97.13 to ?41.73‰ for δD and from ?13.17 to ?4.70‰ for δ¹⁸O. The composition of stable isotopes in the eastern region of this wetland was lower than in the western region. These may have been influenced by uneven evaporation caused by the distribution of aquatic vegetation. During the low-flow season, δD and δ¹⁸O in the more open water with dead aquatic vegetation ranged from ?37.11 to ?11.77‰, and from ?4.25 to ?0.08‰, respectively. This may result from high evaporation rates in this season with the lowest atmospheric humidity. Groundwater fluxes were calculated by mass transfer and isotope mass balance approaches, suggesting that the water sources of the Caohai Wetland were mainly from groundwater in the high-flow season, while the groundwater has a smaller contribution to wetland water during the low-flow season.  相似文献   

Computing derivative information of sequentially coupled subsurface models     

Rafael J. de Moraes  José R. P. Rodrigues  Hadi Hajibeygi  Jan Dirk Jansen 《Computational Geosciences》2018,22(6):1527-1541

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.  相似文献   

Numerical solution and convergence analysis of steam injection in heavy oil reservoirs     

H. Hajinezhad  Ali R. Soheili  Mohammad R. Rasaei  F. Toutounian 《Computational Geosciences》2018,22(6):1433-1444

In this paper, the numerical methods for solving the problem of steam injection in the heavy oil reservoirs are presented. We consider a 3-dimensional model of 3-phase flow, oil, water, and steam, with the effect of 3-phase relative permeability. Interphase mass transfer of water and steam is considered; oil is assumed nonvolatile. We apply the simultaneous solution approach to solve the corresponding nonlinear discretized partial differential equation in the fully implicit form. The convergence of finite difference scheme is proved by the Rosinger theorem. The heuristic Jacobian-Free-Newton-Krylov (HJFNK) method is proposed for solving the system of algebraic equations. The result of this proposed numerical method is well compared with some experimental results. Our numerical results show that the first iteration of the full approximation scheme (FAS) provides a good initial guess for the Newton method. Therefore, we propose a new hybrid-FAS-HJFNK method while there is no steam in the reservoir. The numerical results show that the hybrid-FAS-HJFNK method converges faster than the HJFNK method.  相似文献   

Sensitivity of the Gravity Recovery and Climate Experiment (GRACE) to the complexity of aquifer systems for monitoring of groundwater     

Yashwant B. Katpatal  C. Rishma  Chandan K. Singh 《Hydrogeology Journal》2018,26(3):933-943

The Gravity Recovery and Climate Experiment (GRACE) satellite mission is aimed at assessment of groundwater storage under different terrestrial conditions. The main objective of the presented study is to highlight the significance of aquifer complexity to improve the performance of GRACE in monitoring groundwater. Vidarbha region of Maharashtra, central India, was selected as the study area for analysis, since the region comprises a simple aquifer system in the western region and a complex aquifer system in the eastern region. Groundwater-level-trend analyses of the different aquifer systems and spatial and temporal variation of the terrestrial water storage anomaly were studied to understand the groundwater scenario. GRACE and its field application involve selecting four pixels from the GRACE output with different aquifer systems, where each GRACE pixel encompasses 50–90 monitoring wells. Groundwater storage anomalies (GWSA) are derived for each pixel for the period 2002 to 2015 using the Release 05 (RL05) monthly GRACE gravity models and the Global Land Data Assimilation System (GLDAS) land-surface models (GWSA_GRACE) as well as the actual field data (GWSA_Actual). Correlation analysis between GWSA_GRACE and GWSA_Actual was performed using linear regression. The Pearson and Spearman methods show that the performance of GRACE is good in the region with simple aquifers; however, performance is poorer in the region with multiple aquifer systems. The study highlights the importance of incorporating the sensitivity of GRACE in estimation of groundwater storage in complex aquifer systems in future studies.  相似文献   

The Juchatengo complex: an upper-level ophiolite assemblage of late Paleozoic age in Oaxaca,southern Mexico     

Grajales-Nishimura  José Manuel  Ramos-Arias  Mario Alfredo  Solari  Luigi  Murillo-Muñetón  Gustavo  Centeno-García  Elena  Schaaf  Peter  Torres-Vargas  Ricardo 《International Journal of Earth Sciences》2018,107(3):1005-1031

International Journal of Earth Sciences - The Juchatengo complex (JC) suite is located between the Proterozoic Oaxacan complex to the north and the Xolapa complex to the south, and was amalgamated...  相似文献   

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile   总被引：1，自引：0，他引：1  

Javier Urrutia  Jorge Jódar  Agustín Medina  Christian Herrera  Guillermo Chong  Harry Urqueta  José A. Luque 《Hydrogeology Journal》2018,26(6):1989-2007

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.  相似文献   

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin,northwest China   总被引：1，自引：0，他引：1  

Jingjing Lin  Rui Ma  Yalu Hu  Ziyong Sun  Yanxin Wang  Colin P. R. McCarter 《Hydrogeology Journal》2018,26(5):1559-1572

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?10⁴ m³/year in 2017 to 11,875?×?10⁴ m³/year in 2021, and to 17,039?×?10⁴ m³/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?10⁴ m³/year in 2017 to 1857?×?10⁴ m³/year in 2021, and to 510?×?10⁴ m³/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.  相似文献   

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method     

Tao Chen  Christoph Clauser  Gabriele Marquart  Karen Willbrand  Thomas Hiller 《Hydrogeology Journal》2018,26(6):1903-1916

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.  相似文献   

Tracing submarine groundwater discharge flux in Tolo Harbor,Hong Kong (China)     

Yi Liu  Jiu Jimmy Jiao  Ho Kwan Cheng 《Hydrogeology Journal》2018,26(6):1857-1873

Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m²/d at the study site and 0.3?±?0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25?±?0.07 cm/d) and precipitation (0.45?±?0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m²/d at the study site and 2.42?±?0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12?±?4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed.  相似文献